Smoluchowski's equation for cluster exogenous growth 
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We introduce an extended Smoluchowski equation describing coagulation processes for which 
clusters of mass s grow between collisions with s = As 13 . A physical example, dropwise condensation 
is provided, and its collision kernel K is derived. In the general case, the gelation criterion is 
determined. Exact solutions are found and scaling solutions are investigated. Finally we show how 
these results apply to nucleation of discs on a plane. 



Aggregation processes are of practical interest in many 
fields of science and technology, including chemistry, ma- 
terial sciences, heat transfer engineering, atmosphere sci- 
ences, biology, astronomy, among others [0,0]. In a 
generic aggregation process, clusters of "mass" s, with 
mass distribution N(s,t), can encounter collisions with 
other clusters. Whenever two clusters collide, they can 
stick to form a new cluster with mass conservation. Al- 
though simple numerical models |5]-[5| have proved to be 
a very successful tool in the study of these phenomena, 
much of our theoretical understanding is still to a large 
extent based on Smoluchowski's equation 0, 



d t N(s,t) = \ f* N(s!,t)N(s - sx,t)K(st, s - si) ds x 



-N(s,t) J +oo A^( Sl ,t)K( s , Sl )d Sl , 



(1) 



where the collision kernel K(s±, S2), is the collision rate 
between clusters of mass si and S2, and contains the 
physics of the aggregation process. Smoluchowski's equa- 
tion is obtained in mean-field, neglecting density fluctua- 
tions, and is valid above a model-dependent upper critical 
dimension d c . 

Smoluchowski's equation can be used to study the scal- 
ing properties of a given aggregation model. Dynamic 
scaling corresponds to the fact, observed both in exper- 
iments and in numerical models Q, that the size dis- 
tribution is asymptotically scale-invariant at large time, 
N(s,t) ~ S{ty e f{s/S(t)), where the typical size S(t) 
usually diverges as t z . Smoluchowski's equation can be 
used to determine z and 9, and the profile of the scaling 
function. 

Van Dongcn and Ernst ||[| have extensively studied 
the scaling solutions of Smoluchowski's equation. They 
have shown that collision kernels could be characterized 
by two exponents A and ji, 



K(bx,by) = b x K(x,y), 
K(x,y) ~ (y 



(2) 
(3) 



which determine the small x behavior of f(x). If fi < 0, 
/ is bell-shaped, while for [i > 0, f(x) ~ x~ T . For /1 > 0, 
t = 1 + A, while for fi — 0, r is nontrivial and r < 1 + A. 
In the latter case, such a crucial quantity as the expo- 
nent describing the decay of the total number of clusters 



(n(t) oc t~ z ) is given by z' = z(2 — r) if r > 1, and thus 
depends on r. The determination of r was very difficult 
until a very effective variational method was recently in- 
troduced ||0). 

Despite this reasonably satisfactory state of the art, 
standard Smoluchowski's equation does not describe all 
aggregation processes, and some extensions have been 
provided to include for instance fragmentation of clus- 
ters [0, or source and sink terms ]T^-[T^]. In this Letter, 
we would like to consider another feature, the fact that 
for certain aggregation phenomena, clusters do not grow 
solely due to collisions with other clusters, but also collect 
some mass from the "outside" between collisions events. 
A practical example is dropwise condensation, a daily 
life phenomenon (water condensation on a mirror) which 
has also important consequences in material sciences or 
heat transfer engineering, which motivated most of the 
early studies. Since the seminal work of Beysens and 
Knobler jl5|| , the resulting fascinating droplets patterns, 
or "breath figures", have attracted much interest, and 
simple computer models have been introduced to study 
the kinetics of droplet nucleation fl6|-|l9| , the asymptotic 
surface (or line) coverage [^o|, or the time evolution of 
the "dry" fraction (the surface fraction which has never 
been touched by any droplet) M] (also see the review by 
Meakin fyj). 

In dropwise condensation on a d-dimensional sub- 
strate, individual droplets form D-dimensional hyper- 
spherical caps, and grow by absorption from the vapor 
phase, which leads to a growth law for the droplets ra- 
dius, f oc r w , (uj — —2 for water on a plane Q]), or, 
in terms of the droplet mass s = r D , s = Air, with 
(3 = {D + u — 1)/D. When two droplets overlap, they 
coalesce due to surface tension, to form a new hyper- 
spherical droplet, with mass conservation. Thus, drop- 
wise condensation can be understood as a simple aggrega- 
tion phenomenon, where droplets grow between collisions 
(exogenous growth), as well as through collisions. Such 
a picture is the basis of most simplified computer mod- 
els, such as the droplets growth and coalescence model 
of Family and Meakin |l7j . 

Now, what is the corresponding mean-field kinetic 
equation ? If there were no collisions, the equation would 
simply be a continuity equation expressing the conserva- 
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tion of the number of droplets, 



d t N(s,t)+d s (sN)(s,t) = 0, 



(4) 



with s = s 13 (in the following, we set A = 1). Then, 
the collision rate of two growing droplets of size s\ and 
S2 is, in mean-field, just proportional to the time deriva- 
tive of the cross section er(si, S2) oc (s\ + s^ D ) d , thus 

K(si,s 2 ) oc (s'is 1 / ^ 1 + S2sl /D ^ 1 )(s 1 1 /D + S2 /£ V~\ and 
we obtain the following equation, 



d t N{s,t) + d s ( S PN)(s,t) = 
\ f* N(s u t)N(s - s u t)K( Sl ,s - Sl ,t)d Sl 

- N(s, t) J+°° N(s u t)K(s, Sl ,t)ds u 



with, 



K(x, y) = (x n +yD)(xo + y») 



h\d-l 



(5) 



(6) 



Multiplicative constants were set to 1 by a rescaling of 
t. The expression for K(x,y) is in agreement with an 
early work of Vincent |22], who found a Smoluchowski 
equation for u> — —2, d — 2 and D = 3. However, the 
left hand side of his equation is erroneous since it does 
not conserve the number of particles. The kernel has the 
homogeneity A = (d + uj — 1) / D, as predicted in |l7j] from 
scaling arguments. Eq. (||) for a generic collision ker- 
nel actually describes a general irreversible aggregation 
process with the exogenous growth law s = s' 3 , and we 
would like to study the properties of such an equation. 

Gelation criterion - A very important issue in aggre- 
gation problems is the possible occurrence of gelation. 
Gelation corresponds to the formation of an infinite clus- 
ter at a finite time (in the thermodynamic limit) . Com- 
mon applications are found in food industry for instance. 
Another application in material sciences is the formation 
of fractal aerogels with intriguing physical properties. In 
terms of Smoluchowski's equation, gelation is the phase 
transition associated to the breaking of the mass conser- 
vation through collisions. Adapting standard arguments 
let us consider the net mass flux from clusters with 
s < L towards clusters with s > L, Jr,(t)- From Eq. (||), 
it is easily seen that, 

J L {t)=L 1+ fN(L 1 t) 

+ / dxxN{x,t) / dyK(x,y)N{y,t), (7) 

JO J L-x 

Now, in the absence of an infinite cluster, the mass con- 
servation through collisions requires that Jl (t) — > when 
L — > 00. If gelation occurs at i g , the infinite cluster col- 
lects some mass through exogenous growth and collisions, 
and Ji(t > t g ) has a finite L — > 00 limit, which implies 
that N(s,t) has a slowly vanishing large s tail. Insert- 
ing the ansatz N(s,t) ~ A(t)s~ T for t > t g , we find 
t = max(l + /3, (3 + A)/2). One must have r > 2, since 



the total mass contained in finite clusters must be finite. 
This shows that gelation can occur only if A > 1 or f3 > 1 . 

Exact solutions - For Eq. (|l]), exact solution are known 
only for K = 1, K = x + y, K=xy, but the exis- 
tence of even few exact time-dependent solutions is im- 
portant to check the scaling theories. Here we shall 
provide two exact solutions for Eq. (||), correspond- 
ing to K — 1, with P = and (3=1. To do this, 
we consider the Fourier-Laplace transform of N(s,t), 
Z(z, t) = f + °° e- zs N(s, t) ds. Z(0, t) is the total den- 
sity of clusters n(t). 

For j3 — 0, the Laplace transform of Eq. (||) with 
K = 1 reads, 

d t Z + zZ = Z 2 /2- Z(0,t)Z, (8) 

With n(0) = 1 and Z(z,0) = Z (z), we find Z(0,t) = 
n(t) = 2/(t + 2) and, 



Z(z,t) 



( f + 2 ) 2 (4Z0O) - 3 io JF+W dt ') 



(9) 



Which leads in the scaling regime t — > 00, and for a 
monodispersed initial condition N{s,t) = 5{s — 1) to, 



N(s,t) 



t 2 hit 



e tlnt . 



(10) 



The total mass in the system is, Mi(t) = l + 21n(i+l/2), 
from Mi(t) = -d z Z(0,t). 

For (3=1, the equation for Z is, 



d t Z - zd z Z = Z 2 /2 - Z(0,t)Z. 



(11) 



Once again, n(t) = Z(0, t) = 2/(t + 2), and the solution 
is, 



Z(z,t) 



2Z (ze t ) 



t + 2 (1 - Zo(ze*)) (t + 2) + 2Z (ze t ) ' 



(12) 



For a monodispersed initial distribution, Zq(z) = e~ z , 
Z(z, t) has a pole at zo(t) = — e~* ln(l + 2/t), and we can 
explicitly compute A(s,t), 



N(s,t) 



F expf- S e- t ln(l + ^) ) (13) 



(t + 2) 2 e 



-re te 



which leads in the large time limit to N(s, t) < 

These two exact solutions support the fact that Eq. 
(||) exhibits dynamic scaling, but they also show some 
unusual time dependence of S(t), and a more general 
scaling form, 



AT( s ,t)~y(i)-i/(s/M*)): 



(14) 



where Y(t) does not have the hyperscaling form Y(t) oc 

say. 
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Scaling theory - Now we want to study the scaling solu- 
tions of Eq. (||), with the scaling form of Eq. jl4|), start- 
ing from a monodispersed initial distribution S(s — 1). A 
first remark is that the distribution N(s, t) is at any time 
strictly zero below s (t) = (l + (l-0)t) 1 ^ 1 -^. In phys- 
ical terms, this is due to the fact that the smallest clus- 
ters are the one which have never collided, and thus have 
grown since t = with s = s@ . Consequently, the scaling 
function is strictly zero below x$ = \im(so(t) / S(t)) . xq 
can be zero if S(t) > t 1/(1 - f3 \ but cannot be infinite. 

For kernels with A < 1 and < 1, we can have a 
qualitative understanding of the scaling results to be ex- 
pected. Let us examine Eq. (|j). On the one hand, if 
we switch off the collision term (i.e. the right hand side 
of the equation), the equation describes a set of particles 
which grow in time with s = s 13 , and is associated with 
the size scale S g (t) oc t 1 ^ . On the other hand, if we 
forget the growth term in the left hand side, we are back 
with a standard Smoluchowski equation describing clus- 
tering with mass conservation. The typical size in the 
scaling regime is S c (t) oc fVCi-*), and = 2 §,§. 

Thus, if we switch on both growth and collision, we 
shall observe a "competition" between the dynamic size 
scales corresponding to both processes. If < A, S g (t) ^> 
S c (t), and in the scaling regime we expect S(t) oc S g (t) 
and z = 1/(1 — 0). If > A, on the contrary, the typi- 
cal size of particles increases essentially due to collisions, 
hence S(t) oc S c (t) and z = 1/(1— A). For the marginal 
case A = 0, logarithmic corrections may be observed, and 
are indeed present for the exact solution K = 1, = 0, 

A detailed demonstration of the scaling results is quite 
intricate, since it involves the treatment of many sub- 
cases, and will be published elsewhere EjJ. However the 
main idea is quite simple. Let us define the a-th moment, 
M a (t) = J + °° s a N(s,t)ds. With the scaling form of Eq. 
( |l4| ) , the different terms of Eq. (||) read, 



c/3-l 

d s (s^N)(s : t)^—-(x^fy(x), 



Y 



Collision term 



5 1+A 
Y 2 



(...). 



(15) 



(16) 



(17) 



We look for an asymptotically time independent equation 
consistent with the evolution equation for the total mass, 
obtained by multiplying Eq. (|J) by s and integrating over 
all s, 



Mi = Mp. 



(18) 



We also use the fact that S(t) cannot be much smaller 
than so(t). Results for Y(t) and S(t), are in agreement 
with the qualitative picture above. 



For > A, it is found that S(t) oc so(t), thus z = 
1/(1 - 0), and Y(t) oc S(t) e , with = 2 + A - 0. The 
three terms in the kinetic equation are of the same order 
at large time, and, from the remark above, the scaling 
function is zero below xq = lim(so(t)/S(t)) > 0. For the 
droplets nucleation kernel with d < D, we have (3 > A 
and these results yield z = D/ (1— u>) and = 1 + d/D, in 
agreement with the results for the model of Family and 
Meakin |l7| ]. This scaling theory also shows that f(x) 
cannot diverge at small observed for droplets mod- 
els without nucleation |HJ , and the fact that the scaling 
function is even zero below a finite x$ > is well sup- 
ported by numerical simulations |p3| . 

For < A, S(t) 3> So(t), and the exogenous growth 
term ( |l6| ) is negligible at large time. The scaling equation 
is found to be the same as for standard Smoluchowski's 
equation with the same kernel, 

xf(x) + 2/(1) = f(x) / +o ° f(x x )K(x, x x )d Xl 

~ \ Jo f( x i)f( x - xi)K(xi,x - x\) dxi, (19) 

and the results of (^^] can be applied. The nontrivial 
polydispersity exponents occurring for /i — kernels can 
be computed using the variational method introduced re- 
cently by the present authors p0| . 

For A > 0, we find z = 1/(1 - A), and = 2, the 
total mass in the system being asymptotically conserved 
at large time. 

For the marginal case A = 0, we find different re- 
sults for [L < kernels and /i > kernels. For /i < 
kernels, corresponding to a bell-shaped scaling functions 
(fi < 0), or a nontrivial t < 1 + A (fi = 0), the to- 
tal mass in the system grows logarithmically in the scal- 
ing regime, M\(t) oc hit, we have S(i) oc (t\nt) z : with 
z = 1/(1 - 0), and Y(t) oc S 2 (t)/M 1 (t). For n > 
kernels, with r = 1 + A, there is an additional sublog- 
arithmic correction, leading to M\(t) oc (lnt) ln(lnt), 
S(t) oc (tMi(t)) z , with z unchanged, and Y oc S 2 /M Y 
as before. 

These results are in full agreement with the exact so- 
lution for K = 1 and = 0, both for the asymptotics of 
S,Y, and Mi, and for the scaling function itself. Indeed, 
the scaling function is a pure exponential, just as for the 
exact solution of Eq. ([!]) with K = 1 || . 

Results for systems on the gelling boundary (A = 1 or 
= 1) can also be found and will be reported elsewhere 
p3| . Here, we would like to show how this Smoluchowski 
equation approach can be applied to the nucleation of 
discs on a plane. More generally, we consider the nucle- 
ation model described above, but with d = D. Family 
and Meakin jlTj have simulated growth and coalescence 
of droplets for d = D = 2 and u> = 0.5, and found that 
in contrast with the d < D case, the scaling function is 
polydispersed with a small x divergence f(x) oc x~ T [ p^[ . 
Can we understand this feature ? 
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For this model, the collision kernel is given by Eq. (j|), 
and we have (3 = A + (D — d)/D. Thus, d < D corre- 
sponds to A < (3, while d = D corresponds to A = /3, and 
consequently our scaling theory based on a Smoluchowski 
equation states that the scaling function vanishes below 
a finite xq, for d < D, while / is bell-shaped or polydis- 
persed for d = D, depending on /i. For this kernel, we 
have /i = for lo > 0, and fi = uj/D for u < 0, which 
leads to the cross-over from a bell-shaped scaling func- 
tion for lu < 0, to a small x power law for lo > 0. To 
check this prediction, we performed simulations of the 
droplets growth and coalescence model for d = D = 2 
and various values of lu. These simulations are quite dif- 
ficult since the number of droplets decreases very quickly 
leading to poor statistics and time range limitation. Re- 
sults for the small x behavior of the scaling function for 
different values of lu shown on Fig. [I] are in good agree- 
ment with the expected results. For lu = —3, f(x) has 
a maximum and vanishes faster than any power law as 
x — > 0. When lu — > _ , the position of the maximum 
of the scaling function tends to very quickly, and the 
scaling function crosses over to a power-law at lu = 0. 




-4.0 -3.5 -3.0 
log(s/S(t)) 

FIG. 1. Small x = s/S behavior of the size distributions 
obtained in numerical simulations for different values of the 
growth exponent lu. 

Therefore, the mean-field scaling theory qualitatively 
describes the behavior of the scaling function for d = D, 
and makes for the difference between d < D and d = D. 
Besides, for lo = 0, the mean-field polydispersity expo- 
nent can be computed by the method described in [ p^[ . 
We find r = 1.109 p3[ , which compares well with the 
exponent extracted from the numerics r w 1.2. 

Despite this qualitative success of the Smoluchowski 
equation approach, it should be noticed that the upper 
critical dimension of the model with d = D is proba- 
bly infinite, since we have in mean-field Mi(t) ~ hit, 
whereas Mi(t) is bounded for the actual model, being 
proportional to the surface coverage. 



In conclusion, we have generalized Smoluchowski's 
equation to coagulating systems for which clusters grow 
between collisions with s = s 13 . This equation appears to 
have interesting scaling properties, as found both from 
exact solutions and a general study. This approach re- 
covers scaling results for droplets nucleation models with 
d < D, and qualitatively predicts a transition in the 
shape of the scaling function when varying /3, in agree- 
ment with numerical simulations, despite mean-field lim- 
itations. It would be interesting to see if this kind of 
approach could also be used to describe dropwise con- 
densation with renucleation in empty spaces |jl]|l!|, for 
which nontrivial polydispersity exponents appear. 

We are very grateful to M. Albrecht for a critical read- 
ing of the manuscript. 
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